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Spectral theory and mirror symmetry 

Marcos Marino 


Abstract. Recent developments in string theory have revealed a surprising 
connection between spectral theory and local mirror symmetry: it has been 
found that the quantization of mirror curves to toric Calabi-Yau threefolds 
leads to trace class operators, whose spectral properties are conjecturally en¬ 
coded in the enumerative geometry of the Calabi—Yau. This leads to a new, 
infinite family of solvable spectral problems: the Fredholm determinants of 
these operators can be found explicitly in terms of Gromov-Witten invari¬ 
ants and their refinements; their spectrum is encoded in exact quantization 
conditions, and turns out to be determined by the vanishing of a quantum 
theta function. Conversely, the spectral theory of these operators provides a 
non-perturbative definition of topological string theory on toric Calabi—Yau 
threefolds. In particular, their integral kernels lead to matrix integral rep¬ 
resentations of the topological string partition function, which explain some 
number-theoretic properties of the periods. In this paper we give a peda¬ 
gogical overview of these developments with a focus on their mathematical 
implications. 


1. Introduction 

Mirror symmetry has played a fundamental role in the interface of theoretical 
physics and mathematics, and has led to beautiful developments in string theory 
and enumerative geometry. Although mirror symmetry was originally formulated 
for compact Calabi-Yau (CY) manifolds, it can be extended to the so-called local 
case [KKV, CKYZ], which involves toric, hence non-compact CY manifolds. Lo¬ 
cal mirror symmetry is in a sense even richer than its compact counterpart, since it 
is related to many other fields of mathematical physics, like Chern-Simons theory, 
integrable systems, supersymmetric gauge theory, and random matrix theory. Very 
often, these connections have been obtained by studying topological string theory, 
the physical theory underlying mirror symmetry. In some circumstances, this theory 
can be described by other, “dual” theories, which look in principle very different, 
and this physical equivalence leads to a non-trivial mathematical equivalence. For 
example, the large N duality obtained by Gopakumar and Vafa in [GVl] relates 
topological string theory on the resolved conifold to Chern-Simons theory on the 
three-sphere, and this equivalence eventually led to the theory of the topological 
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vertex [AKMV], which gives a complete algorithm to obtain the Gromov-Witten 
invariants of toric CY manifolds in terms of Chern-Simons knot invariants. 

The spectral theory of self-adjoint operators plays a more established role in 
physics, due to its connection to Quantum Mechanics. In recent years, there have 
been many interesting developments in the arena of exactly solvable spectral prob¬ 
lems. For example, in some one-dimensional models, exact quantization conditions 
incorporating non-perturbative effects at all orders have been found (see for exam¬ 
ple [ZJJ, Vo3, Al]). Some of these solvable models are closely related to integrable 
systems [DDT] and two-dimensional CFT [BLZ]. 

The purpose of this review paper is to give a summary of some recent results 
in string theory which suggest a deep connection between these two separate sub¬ 
jects. There is now growing evidence that, given a toric CY manifold, one can 
associate to it trace class operators whose spectral properties are encoded in the 
enumerative invariants of the underlying manifold. This correspondence is quite 
precise and leads to conjectural, explicit formulae for the Fredholm determinants of 
these operators, from which one can derive exact quantization conditions for their 
spectrum. From the point of view of spectral theory, this connection provides a new 
family of solvable operators with beautiful properties. From the point of view of 
enumerative geometry, it gives a new meaning to the generating functionals encod¬ 
ing the enumerative invariants of local CY geometries: they appear as asymptotic 
expansions of quantities defined by the spectral theory of the corresponding quan¬ 
tum operators. For example, one can obtain the genus expansion of the topological 
string free energy as an asymptotic expansion of the “fermionic” spectral traces of 
these operators. 

The correspondence between spectral theory and topological strings which we 
will describe can be regarded in many ways as a large N string/gauge duality, 
like those appearing in the AdS/CFT correspondence. It relates a simple quantum- 
mechanical problem on the real line to a string propagating on a toric CY manifold. 
The weak coupling regime of the quantum problem, when the Planck constant is 
small, corresponds to the regime in which the string coupling constant is strong. 
The genus expansion of the topological string emerges in a’t Hooft-like limit of the 
quantum mechanical problem, which can be regarded as a rigorous non-perturbative 
definition of topological string theory on these backgrounds. 

This review paper is organized as follows. In section 2, we present a simple 
example of the general type of spectral problems we will deal with, and we present 
both the approximate solution obtained with elementary methods and the exact 
solution conjectured in the recent literature. In section 3 we present the general 
conjecture formulated in [GHM] connecting spectral theory and topological string 
theory on toric CY manifolds. The emphasis is on how string theory solves the 
spectral problem. In section 4 we explain how, reciprocally, the spectral problem 
leads to new insights on topological string theory, providing in particular a non- 
perturbative definition of the genus expansion. Finally, in section 5, we conclude 
with some interesting problems for the future. 
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2. A problem in spectral theory 

2.1. A quantum curve. Let us start by formulating a sharp question con¬ 
cerning the spectrum of an interesting operator on the real line. This will hopefully 
provide a motivation for the type of problems that we would like to address in this 
overview paper. 

Let X, y be Heisenberg operators on the real line, satisfying the commutation 
relation 

(1) [x,y] = i;i. 

Let us consider the following operator, 

(2) 0p2 = e=<-f e^-f e-^'-L 

Why look at this beast? The reason is that it can be obtained by “quantizing” the 
mirror curve to the CY threefold known as “local P^”, which is the total space of 
the canonical bundle of P^: 

(3) A = 0(-3)^P^ 

This is probably the most studied example among all non-compact CY manifolds. 
As we will briefly review later on, the mirror of this manifold is encoded in an 
elliptic curve in (D* x (D*, which can be written in the form 

(4) + K = t). 

Here, n parametrizes the moduli space of complex structures of the curve. We 
can now perform the standard Weyl quantization of the first three terms of the 
curve (4), i.e. we promote x and y to Heisenberg operators and we use the Weyl 
ordering prescription. This gives the operator (2). It involves the exponentiated 
“coordinate” and “momentum” operators 

(5) X = e^ Y = e^, 
which are self-adjoint and they satisfy the Weyl algebra 

(6) XY = gYX, 
where 

(7) q = 

The domain of the operator X, D{X), consists of functions ^ L^(R) such that 

(8) e^i;{x) e l2(R). 

Similarly, the domain of Y, D{Y), consists of functions ipix) G L^(R) such that 

(9) ey^(y) G l2(R), 

where 

(10) ^^(y) = J 
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n 

En 

0 

2.56264206862381937 

1 

3.91821318829983977 

2 

4.91178982376733606 

3 

5.73573703542155946 

4 

6.45535922844299896 


Table 1. Numerical spectrum of the operator (2) for n = 
0,1, • • • ,4, and h = 27r. 


is the wavefunction in the y representation, which is essentially given by a Fourier 
transform. The condition (9) can be translated into a condition on requiring 

analyticity on strips. The domain of the operator 0p2 should be chosen appropri¬ 
ately, by considering functions such that 0 p 2 i/! G L^(R). 

As we will review in the next section, it was shown in [KM] that the operator 
0p2 has an inverse which is positive-definite and of trace class in L^(R), so its 
spectrum (which depends on h) consists of an infinite set of positive eigenvalues 
e^", n = 0,1,---. This spectrum can be determined numerically with standard 
methods, as first pointed out in [HW]. First, one should choose an appropriate 
orthonormal basis of functions - of L^(R), in the domain of 0p2 (for 

example, the eigenfunctions of the harmonic oscillator will do). Then, the e^" are 
the eigenvalues of the infinite-dimensional matrix 

( 11 ) Mij = f,j = 0 , 1 ,---. 

They can be obtained numerically by truncating the matrix to very large sizes. 
Boundary effects can be partially eliminated with standard extrapolation methods, 
and one can find very accurate values for the En. For h = 27r, the results for the 
very first eigenvalues are listed in Table 1. We can now ask the following question: 
can these eigenvalues be determined analytically! Usually, this is hopeless, since 
there are very few operators whose spectral properties can be determined exactly. 
However, a conjecture put forward in [GHM] affirms that in this case it can be 
done: there is an exact, conjectural expression for the Fredholm determinant of the 
inverse operator 

(12) PP2=Op2, 

from which one can deduce an exact quantization condition for the spectrum^. 
Before presenting this conjecture, let us see what can be said about the spectrum 
of the above operator with elementary methods. 

2.2. A first approach to the spectrum. Although the above spectral prob¬ 
lem is not of the standard form found in Quantum Mechanics, it is possible to use 
conventional approximation techniques to gain some insight. In particular, one can 
try to use the WKB approximation. Let us briefly review this method. Let us 
assume that we have a classical Hamiltonian of the form, 

(13) H{x,y) = ^+V{x), 

^The exact quantization condition for this particular operator and h = 2 tt was first written 
down in [HW] by generalizing a result of [KaMa] for a different CY geometry. 
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where y is interpreted as the momentum, and V(x) is a potential supporting bound 
states. We will assume for simplicity that V{x) is a confining potential, i.e. that 
V{x) —)■ +00 for |a;| —> oo. After quantization, one hnds a quantum Hamiltonian 
H with a discrete spectrum of energies En, n = 0,1, • • •. The Bohr-Sommerfeld 
quantization condition gives an approximate quantization condition which deter¬ 
mines the spectrum in the limit of small h or large quantum numbers n ^ 1. It 
can be formulated geometrically as follows. One considers the plane curve 

(14) H{x,y)=E, 

where E is interpreted as the energy of the system, and a differential on this curve 
given by 

(15) A = y{x)dx. 

In simple situations, there is a B-cycle in the curve (14) which corresponds to 
periodic motion in the potential, between two turning points determined by the 
equation of the curve. The integral of A over this cycle gives the volume of the 
region TZ{E) in phase space with energy less or equal than E, 

(16) 7^(A) = {(x, 2 /) e : H{x, y) < E}. 

We will denote this volume by volo(i?). The Bohr-Sommerfeld quantization condi¬ 
tion reads 

(17) volo(i?) = X = 2Trh + 2^ ’ n = 0,1, • • • . 

It can be also interpreted as saying that each cell of volume 27r/i in Tl{E) corresponds 
roughly to a quantum state. 

It turns out that, although our quantum operator (2) is not of the standard form 
(13), we can still use the Bohr-Sommerfeld condition. The classical counterpart 
of the operator 0 p 2 is the function on the phase space given by the first three 
terms in (4), with the symplectic form w = dx A dy, 

(18) Op^{x,y) =e^ +ey +e-^-y. 

The analogue of the hyperelliptic curve (14) is now, 

(19) C>p2(x,2/) = e®, 

and the function (18) defines a compact region in phase space 

(20) n{E) = {{x,y) e R2 : Op 2 (x,y) < e^}, 

which is the analogue of (16). Note that we are implicitly interpreting Op 2 as the 
exponential of a classical Hamiltonian. The Bohr-Sommerfeld quantization con¬ 
dition for this problem is again given by (17). The period of A can be computed 
explicitly, but since the resulting quantization condition is only valid at large ener¬ 
gies, it is worth to simplify it, as follows. At large E, the region TZ{E) has a natural 
mathematical interpretation: it is the region enclosed by the tropical limit [BIMS] 
of the curve 

( 21 ) + ey + e-^-y - = 0, 

and limited by the lines 

(22) X = A, 


y = E, 


X + y + E = 0, 
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see Fig. 1 for an illustration when E = 15. In this polygonal limit, we find 

(23) volo(i?) « 

and we obtain from (17) the following approximate behavior for the eigenvalues, at 
large quantum numbers, 

(24) n» 1. 

O 

It can be seen, by direct examination of the numerical spectrum, that this rough 
WKB estimate gives a good approximation to the eigenvalues of the operator (2) 
when n is large. The estimate (24) has been rigorously proved in [LST]. 



Figure 1. The region TZ{E) for E = 15. 


2.3. An exact quantization condition. An important question is whether 
the Bohr-Sommerfeld quantization condition obtained above can be improved. One 
possible strategy is to consider the WKB quantization scheme to all orders in h. 
The all-orders WKB quantization condition was derived for standard Schrddinger 
operators by Dunham in [D]. It can be adapted to operators of the form (12) 
by using a general formalism due to Voros [Vol], or by using the extension of 
the WKB method to difference equations in [DM]. For the operator (2), this all¬ 
orders WKB method was studied in [ACDKV]. As we will explain in detail in 
the next section, the main (conjectural) result of [ACDKV] is that the resulting 
quantization condition can be expressed in terms of enumerative invariants of the 
CY (the so-called NS limit of the refined topological string.) However, as pointed 
out in [KaMa], this can not be the full answer to the spectral problem, since it leads 
to unacceptable poles for values of h of the form h = qir, where q G Q. Physically, 
one has to take into account the contribution from quantum-mechanical instantons, 
which are known to correct the all-orders WKB quantization condition [BPV]. 

It turns out that in the case h = 27r, the exact quantization condition is re¬ 
markably simple. Let us consider the Picard-Fuchs equation, 

(25) (6»3-3z(36»-H2)(36»-bl)6i)n = 0, 

where 


d 
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This equation has three independent solutions at z = 0: a trivial, constant solu¬ 
tion, a logarithmic solution Wi{z) and a double logarithmic solution tx72(z). If we 
introduce the power series. 


(27) 


= 2.3——(-z)^ 

j>l 


^2{z) = XI {^^(3^') - i’U + 1)} 


j> 


t j! r(i + jT 


where is the digamma function, we have that 


(28) 


Wi{z) = log(z) -h t77i(z), 

W 2 iz) = log^(z) -I- 2^1 (z) log(z) -I- W 2 {z). 


It is easy to see that the series in (27) have a radius of convergence |z| = 1/27. The 
solutions 'Ci7i,2(z:) can be written in terms of hypergeometric and Meijer functions. 
Let us now introduce the following function, 

(29) m = 1 

oTT LU ^ '2' ^ 

where 

(30) z = e-3^. 

Then, the conjecture of [GHM] implies that the spectrum of (2), for h = 2 tt, is 
determined by the condition 

(31) ^{E)-^=n+^, n = 0,1,---. 

It can be easily seen that the values of En determined in this way agree with the 
ones in Table 1. A graphical representation of the function ^{E) — 1/4, showing the 
different energy levels, can be found in Fig. 2. 



Figure 2. The function ^{E) — 1/4, as a function of E, determin¬ 
ing the energy spectrum of the operator (2) for h = 27r. 

What is the content of the quantization condition (31)? The function ^{E) has 
a convergent power series expansion of the form, 

Qp2 °° oo 

^E) = ^+eJ2 0-6-3"^ + X 

n—1 n—1 


( 32 ) 











MARCOS MARINO 


(33) 


m- 


When E is large, we can approximate 

^ 9E^ 

4 Stt^ ’ 

and we recover the Bohr-Sommerfeld estimate (24) of the previous section. The 
constant 1/4 is a correction which can be easily obtained with the next-to-leading 
WKB method [KaMa, HW, GHM], and then we have an inhnite series of expo¬ 
nentially small corrections in the energy. As we will explain later on, they have a 
mixed physical origin: they combine higher order, perturbative WKB corrections 
with non-perturbative, instanton-type corrections. 

The reader familiar with local mirror symmetry has surely recognized the equa¬ 
tion (25): it is the Picard-Fuchs equation governing the periods (28) of the mirror 
of local P^, which fully determine the genus zero Gromov-Witten invariants of this 
CY. The variable z is a modulus for the mirror curve, and it is related to k in (4) 
by 

(34) z = K 


-3 


Therefore, the corrections to the Bohr-Sommerfeld quantization condition are gov¬ 
erned by the topological string on local P^, whose mirror is precisely the algebraic 
curve (4)! As we will see in the next section, this is not a coincidence, but a general 
story which seems to hold for all toric CY threefolds: the spectral theory of the 
operators obtained by quantizing their mirror curves is solved by their enumerative 
invariants. 


3. Prom topological strings to spectral theory 

3.1. Mirror symmetry and topological strings. In this section we will 
make a quick summary of the relevant background of mirror symmetry and topo¬ 
logical string theory. Useful reviews of these subjects might be found in [CK, H+]. 

Mirror symmetry relates pairs of CY threefolds, X, X, in such a way that the 
enumerative geometry of X is reformulated in terms of the deformation of complex 
structures of X. There are various types of enumerative invariants of X. Let us 
consider holomorphic maps from a Riemann surface of genus g to the CY X, 

(35) f:Eg^X. 

Let [5^] G H2{X,Wi)^ i = I,-- - ,s, be a basis for the two-homology of X, with 
s = b 2 {X). The maps (35) are classified topologically by the homology class 

S 

(36) /*[(S,)]=E dip,] G H2{X,I.), 

i=l 

where di are integers called the degrees of the map. We will put them together 
in a degree vector d = (di,--- ,ds). The Gromov-Witten invariant at genus g 
and degree d, which we will denote by N^, “counts” (in an appropriate way) the 
number of holomorphic maps of degree d from a Riemann surface of genus g to 
the CY X. Due to the nature of the moduli space of maps, these invariants are in 
general rational, rather than integer, numbers; see for example [CK] for rigorous 
definitions and examples. 

The Gromov-Witten invariants at fixed genus g but at all degrees can be put 
together in generating functionals Fg{t), usually called genus g free energies. These 
are formal power series in e“*L * = Ir ’ ’ ? S) where ti are the Kdhler parameters 
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of X. More precisely, the ti are flat coordinates on the moduli space of Kahler 
structures of X. It is convenient to add to these generating functionals polynomial 
terms which appear naturally in the study of mirror symmetry and topological 
strings. In this way, we have, at genus zero, 

1 ® 

(37) Fo{t) — - Qijktitjtk NqB 

d 

In the case of a compact CY threefold, the numbers aijk are interpreted as triple 
intersection numbers of two-classes in X. At genus one, one has 

S 

(38) = + 

i=l d 

In the compact case, the coefficients bi are related to the second Chern class of the 
CY manifold [BCOVl], At higher genus one finds 

(39) Fg(t) = Cg+5]<e-‘i‘, g>2, 

d 

where Cg is a constant, called the constant map contribution to the free energy 
[BCOV2]. It turns out that these functionals have a physical interpretation as the 
free energies at genus g of the so-called type A topological string on X. Roughly 
speaking, this free energy can be computed by considering the path integral of 
a string theory on Riemann surfaces or “worldsheets” of genus g. This makes it 
possible to use a large amount of methods and ideas of physics in order to shed 
light on these quantities. 

Although the above generating functionals are in principle formal generating 
functions, they have a common region of convergence near the large radius point 
ti —>■ oo. The total free energy of the topological string is formally defined as the 
sum, 

(40) F^^ (t,g,) = ^ 

g>0 g>0 d 

where 

s s 

(41) F‘'P\t,gs) = -^ hU + Y 

The superscript “WS” refers to worldsheet instantons, which are counted by this 
generating functional. The variable gs, called the topological string coupling con¬ 
stant, is in principle a formal variable keeping track of the genus. However, in string 
theory this constant has a physical meaning, and measures the strength of the string 
interaction. When g^ is very small, only Riemann surfaces of low genus contribute 
to a given quantum observable. On the contrary, if gs is large, the contribution of 
higher genus Riemann surfaces becomes very important. 

The convergence properties of the formal series (40) are less understood, but 
there is strong evidence that, for fixed U in the common region of convergence, 
the numerical series Fg{t) diverges factorially, as (2^)! (see [M2] and references 
therein.) Therefore, the total free energy (40) does not define in principle a function 
of gs and t^. In a remarkable paper [GV2], Gopakumar and Vafa pointed out that 
the series (40) can be however resummed order by order in exp(—t,), at all orders 
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in gs- This resummation involves a new set of enumerative invariants, the so- 
called Gopakumar-Vafa invariants n'^. Out of these invariants, one constructs the 
generating series 



g>0 d w—1 

and one has, as an equality of formal series, 


(43) (t, g^) = (t, g,) + (t, g ^). 


The Gopakumar-Vafa invariants turn out to be integers, in contrast to the original 
Gromov-Witten invariants. One can obtain one set of invariants from the other by 
comparing (40) to (43), but there exist direct mathematical constructions of the 
Gopakumar-Vafa invariants as well, see [PT] . 

The mirror manifold X to X has s complex deformation parameters Zi, i = 
1, • • • , s, which are related to the Kahler parameters of X by the so-called mirror 
map ti{z). At genus zero, the theory of deformation of complex structures of X 
can be formulated as a theory of periods for the holomorphic (3, 0) form of the CY 
X, which vary with the complex structure moduli Zi. As first noted in [CGPO], 
this theory is remarkably simple. In particular, one can determine a function Fq(z) 
depending on the complex moduli, called in this context the prepotential, which is 
the generating functional of genus zero Gromov-Witten invariants (37), once the 
mirror map is used. The theory of deformation of complex structures also has a 
physical realization in terms of the so-called type B topological string. However, 
for general GY manifolds, a precise definition of this theory at higher genus is 
still lacking. In practice, one often uses the description of the B-model provided 
in [BCOV2]. This provides a set of constraints for a series of functions Fg{z), 
5 > I, known as holomorphic anomaly equations. After using the mirror map, 
these functions become the generating functionals (39). 

One of the consequences of mirror symmetry and the B-model is the existence 
of many possible descriptions of the topological string, related by symplectic trans¬ 
formations. The existence of these descriptions can be regarded as a generalization 
of electric-magnetic duality. In this way, one finds different “frames” in which the 
topological string amplitudes can be expressed. Although these frames are in prin¬ 
ciple equivalent, some of them might be more convenient, depending on the region 
of moduli space we are looking at. The topological string amplitudes Fg(t) written 
down above, in terms of Gromov-Witten invariants, correspond to the so-called 
large radius frame in the B-model, and they are appropriate for the so-called large 
radius limit Re(tj) ^ 1. The topological string free energies in different frames are 
related by a formal Fourier-Laplace transform, as explained in [ABK]. 

Although mirror symmetry was formulated originally for compact GY three¬ 
folds, one can extend it to the so-called “local” case [KKV, CKYZ]. In local 
mirror symmetry, the CY X is taken to be a toric CY manifold, which is necessar¬ 
ily non-compact. The theory of deformation of complex structures of the mirror X 
is encoded in an algebraic curve of the form 


(44) 


IV(e“^,e^) = 0, 


where the variables appear naturally exponentiated. Local mirror symmetry is 
considerably simpler than full-fledged mirror symmetry. On the A-model side, the 
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enumerative invariants can be computed algorithmically in various ways, either by 
localization [Ko, GP, KZ], or by using the so-called topological vertex [AKMV], 
At the same time, the theory of deformation of complex structures of X can be 
simplified very much. At genus zero, one should consider the periods of the differ¬ 
ential 


(45) A = y{x)dx 

on the curve (44) [KKV, CKYZ]. The mirror map and the genus zero free energy 
To(t) in the large radius frame are determined by making an appropriate choice of 
cycles on the curve, Oi, /3i, z = 1, • • • , s, and one finds 

(46) ^ * = lc",s. 

In general, s > gs, where gs is the genus of the mirror curve. Also, in the local case, 
the type B topological string can be formulated in a more precise way, by using the 
topological recursion of [EOl], in terms of periods and residues of meromorphic 
forms on the curve (44) [Ml, BKMP]. In particular, mirror symmetry can be 
proved to all genera [E02, FLZ] by comparing the definition of the B model in 
[Ml, BKMP] with localization computations in the A model. 

Example 3.1. Local P^. In the case of the local CY given in (3), the genus 
zero free energy can be obtained by mirror symmetry as follows. The complex 
deformation parameter in the mirror curve is the parameter z appearing in (4). 
The mirror map is given by 


(47) 


t = -wi{z), 


where tui(z) is the power series in (28). Then, Fo(t) is defined, up to a constant, 
by 


(48) 


dFp _ W2(z) 

dt ~ 6 


where '072 (z) is the other period in (28). One then finds, after fixing the integration 
constant appropriately, 


(49) 




18 
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From here one can read the very first Gromov-Witten invariants, like for example 
^ = 3, see [CKYZ] for more details. □ 


Another interesting feature of the local case is that it is possible to define a more 
general set of enumerative invariants, and consequently a more general topological 
string theory, known sometimes as the refined topological string. This refinement 
has its roots in the instanton partition functions of Nekrasov for supersymmetric 
gauge theories [N] . Different aspects of the refinement have been worked out in for 
example [IKV, KW, HK]. In particular, one can generalize the Gopakumar-Vafa 
invariants to the so-called refined BPS invariants. Precise mathematical definitions 
can be found in [CKK, NO]. These invariants, which are also integers, depend on 
the degrees d and on two non-negative half-integers, Jl, Jr, or “spins”. We will de¬ 
note them by . The Gopakumar-Vafa invariants are particular combinations 
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of these refined BPS invariants, and one has the following relationship, 

(50) X] ^ , 

iijH 9>0 


where g is a formal variable and 

(51) X]{q) 



is the SU(2) character for the spin j. We note that the sums in (50) are well- 
defined, since for given degrees d only a finite number of jji, g give a non-zero 
contribution. 

Out of these refined BPS invariants, one can define the Nekrasov-Shatashvili 
(NS) free energy. 


^ S S 

^) = ^ X! O-ijktitjtk + X] 


(52) 




2=1 


+ 

3 l, 3 r 


3 R 


sinf^(2j^ + l)sin'f(2jfi + l) ^_^^.^ 
2w2 sin^ ^ ® 


In this equation, the coefficients aijk are the same ones that appear in (37), while 
can be obtained by using mirror symmetry as in [KW, HK], The free energy 
(52) is not the most general generating functional for the refined BPS invariants, 
but involves a particular combination thereof, which defines the so-called NS limit 
of the refined topological string. This limit was first discussed in the context of 
gauge theory in [NS]. By expanding (52) in powers of H, we find the NS free 
energies at order n, F^^{t), as 


(53) F^^{t,h) = X]FNS(t)fi2»-i. 

n=0 

The expression (52) can be regarded as a Gopakumar-Vafa-like resummation of 
the series in (53). An important observation is that the first term in this series, 
is equal to Tb(t), the standard genus zero free energy. Note that the term 
involving the coefficients contributes to F;(^®(t). 


3.2. Quantizing curves. The semiclassical limit of many interesting prob¬ 
lems in quantum physics turns out to be described by a plane curve, together with 
a choice of a meromorphic differential on it. Perhaps the simplest example is the 
Bohr-Sommerfeld quantization condition (17) for the one-dimensional Schrodinger 
equation. A more complicated example is the planar limit of random matrix mod¬ 
els, which is typically encoded in a spectral curve which captures the information 
on the limiting distribution of eigenvalues (here, the semiclassical limit is the large 
N limit). Another example which has received much attention recently is the semi¬ 
classical limit of SL(2,(D) Chern-Simons theory on the complement of a knot. In 
this case, the relevant algebraic curve is the A-polynomial of a knot [Gj. 

The genus zero free energy of topological string theory can be formally regarded 
as a semiclassical limit, in which the role of h is played by the string coupling 
constant in the limit —>■ 0, the dominant term in the genus expansion of (40) 
is indeed the genus zero free energy. As we explained in (46), this semiclassical 
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limit can be obtained by computing periods of the meromorphic differential (15) 
on the mirror curve. 

It is tempting to think that, in a situation where an algebraic curve encodes 
the classical limit of a quantum problem, the quantum corrections are obtained by 
considering a “quantum” version of the curve. This has led to a rich literature on 
“quantum curves,” where the original curve is typically promoted to a differential 
operator (see for example [No, GS] for recent reviews). In the case of topolog¬ 
ical string theory, an approach based on quantum curves was first proposed in 
[ADKMV], where it was suggested that the higher genus free energies of topo¬ 
logical string theory (which can be regarded as “quantum corrections”) might be 
obtained by quantizing the mirror curve in an appropriate way. Building on work 
on supersymmetric gauge theories and quantum integrable systems [NS], it was 
later shown in [ACDKV] that mirror curves can be formally quantized by using 
the WKB approximation (this was first pointed out in [MM], in the context of 
Seiberg-Witten curves). However, the quantum corrections obtained in this way 
do not give the conventional higher genus Gromov-Witten invariants, but rather 
the NS free energy (53). 

It is well-known that the quantization of a classical system is in general not 
unique: there can be different quantizations of the same function, leading to the 
same classical limit but differing when H 0. The quantization of an algebraic 
curve is also plagued with ambiguities, for at least three different reasons. First of 
all, the quantization of a curve inherits the ordering ambiguities of any quantization 
procedure. Second, algebraic curves are naturally defined in the complex domain, 
while their quantum versions typically involve a choice of reality conditions which 
can be made in different ways. Finally, when an algebraic curve is promoted to an 
operator, one has to specify the Hilbert space of wavefunctions where this operator 
acts, and in particular the boundary conditions satisfied by such wavefunctions. In 
most of the recent literature on “quantum curves,” these issues are not addressed 
explicitly, and quantum curves are only studied at the level of formal WKB expan¬ 
sions, as in [ACDKV]. This perturbative approach has led to many interesting 
results, but does not give a non-perturbative formulation of the quantum problem. 

In [KaMa] it was pointed out that the mirror curves of toric CY threefolds can 
be quantized in such a way that the resulting operators have a discrete spectrum 
(we have already seen an example in section 2, in the case of local P^.) It was 
then conjectured in [GHM] that the quantization of mirror curves leads to posi¬ 
tive definite, trace class operators on L^(R). This was proved in [KM] in many 
examples. In the scheme proposed in [GHM, KM], the ambiguities arising in the 
quantization of algebraic curves are solved in a natural way: ordering ambiguities 
are fixed by using Weyl quantization, which is particularly well suited for exponen¬ 
tiated position and momentum operators. Reality conditions are chosen in such a 
way that the classical regions (20) in phase space are compact. Finally, the Hilbert 
space where the operators act is simply L^(R). 

In the rest of this paper, we will focus for simplicity on toric (almost) del Pezzo 
CY threefolds. These CYs are defined as the total space of the canonical bundle 
on a toric (almost) del Pezzo surface S, 

(54) A = 0{Ks) S. 

They are sometimes called “local S.” For example, if S' = P^, the total space of 
its canonical bundle will be called local P^. Examples of toric del Pezzos include. 
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besides P^, the Hirzebruch surfaces F„, n — 0,1,2, and the blowups of at n 
points, denoted by S„, for n = 1, 2,3 (note that Fi = Bi, and that Fq = P^ x P^). 
The main simplifying characteristic of these manifolds is that their mirror curve has 
genus one. This makes their analysis much simpler. We hasten to add that there 
is a very interesting generalization of all the considerations in this review paper to 
mirror curves of higher genus [CGM2]. In this case, a mirror curve of genus gs 
leads in general to gs trace class operators. 

By standard results in toric geometry (see for example [HKP, CR]), toric, 
almost del Pezzo surfaces can be classified by reflexive polyhedra in two dimensions. 
The polyhedron Ag associated to a surface S is the convex hull of a set of two- 
dimensional vectors 

(55) !/(*) = i = 1, • • • ,s-f 2, 


together with the origin, see Fig. 3 for the example of local P^. In order to construct 
the total space of the canonical bundle over S, we have to consider the extended 
vectors 


(56) 


= (1,0,0), 

I7W = (l,iz«,i.W), 


They satisfy the relations 


i — I, • • • , s -|- 2. 


(57) 


s+2 


2 = 0 


= 0 , 


ex 1, * * * , s. 


where Qf is a matrix of integers (called the charge matrix) which characterizes the 
geometry. 

The construction of the mirror geometry to (54) goes back to Batyrev, and 
it has been recently reviewed in [CR], to which we refer for further details. In 
order to write down the equation for the mirror curve of (54), we note that the s 
complex parameters of the mirror can be divided in two types: one “true” modulus 
K (for a curve of genus one) and a set of “mass” parameters i = I,-- - ,s — I 
[HKP, HKRSj. In terms of these variables, the mirror curve for a local del Pezzo 
CY threefold can be written as. 


(58) 
where 

(59) 


W{e\ey)=Os{x,y) + fe = 0, 


S+2 

Os{x,y) = ^exp 


ul X + 1x2 y + 


m)) 


and fi{^) are suitable functions of the parameters 


Example 3.2. The simplest case of a local del Pezzo CY is local P^. In this 
case, we have s = 1. The vectors (55) are given by 

(60) = (1,0), +) = (0,I), +) = (-!,-!). 


We show these vectors, together with the origin and their convex hull Ap 2 , in Fig. 3. 
In this geometry there is one complex deformation parameter x, and the function 
Op 2 {x,y) is given by 

(61) 


C>p 2 {x,y) = + ey + e-^-y 
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Figure 3. The vectors (60) defining the local geometry, to¬ 
gether with the polyhedron Ap 2 (in thick lines) and the dual poly¬ 
hedron (in dashed lines). 


s Osix,y) 

P^ +ey + 

Fq 

+ey + Q-y 

Fi 

+ey + e-^-y + 

F2 

+ey + e-'^^-y + ^e'^ 

B2 

+ey + e-^-y + Cie"*' -k ^26"“ 

B3 

+ey + e-^-y + ^le-^ + ^26"^ -k ^se^+y 


Table 2. The functions Osix,y) associated to some local del 
Pezzo CYs. 


This is the curve (4). Note that, as shown in Fig. 3, the dual polyhedron to Ap 2 is 
precisely the region in Fig. 1. □ 

Example 3.3. The previous example can be generalized by considering the 
canonical bundle over the weighted projective space P(l, m,n), where m,n S S>o- 
This is not a smooth manifold, but it can be analyzed by using extensions of 
Gromov-Witten theory, see for example [BC] for a study of the case n = 1. The 
vectors are in this case 

(62) p^^^ = (l,0), u^^^ = (0,1), = (—m, —n), 

and the function appearing in the mirror curve (58) is given by 

(63) {x,y)=e^+ey + . 

Some of these geometries can arise as degeneration limits of toric del Pezzos. For 
example, the mirror curve to local F 2 is characterized by the function 

(64) Op, (x,y)=e"+e^ + e-2“-^+^e-", 

and when ^ = 0 we recover the geometry (63) with m = 2 and n = 1. □ 

Some further examples of functions obtained by quantization of local del Pezzos 
can be found in Table 2. Details on the corresponding geometries can be found in 
for example [HKP]. 
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The “quantization” of the mirror curve (58), in the case of local del Pezzo 
CYs, is based on the promotion of the function Osix,y) to an operator, which 
will be denoted by Os- First, we promote a;, y to self-adjoint Heisenberg operators 
X, y on the real line, satisfying the commutation relation (1). This involves a 
choice of reality conditions for the complex variables x, y. Then, we apply Weyl’s 
quantization to Osix, y). In particular, this fixes possible ordering ambiguities. We 
recall that, in Weyl’s quantization, we have 

(65) ^ e™+"y. 

As noted in [GHM], instead of focusing on the operator Os (which is not of trace 
class), one should rather consider its inverse 

( 66 ) PS = 

acting on L^(R). It was conjectured in [GHM] that the operators ps are of trace 
class and positive definite, provided appropriate positivity conditions are imposed 
on the mass parameters appearing in the mirror curves. 

In order to motivate this assertion, note that, if e'®’* are the eigenvalues of Os, 
the eigenvalues of ps are e“®". It is easy to extend the Bohr-Sommerfeld estimate 
(24) to all the operators Os (see [GHM]): we can obviously define an available 
region in phase space, TZ{E), similar to what we did in (20), but involving this time 
the classical curve Os- The volume of this region behaves, at large E, as 

(67) vo1o(F;) « C'F;^ A>1. 


The coefficient C can be easily computed by considering the tropical limit of the 
classical curve, where the region TZ{E) becomes a polygon (in fact, it is the dual 
polyhedron to As). We then find the estimate 

(68) En « n > 1. 

This heuristic argument indicates that the spectral traces of the operators ps, 

(69) Trp| = ^e-'^-, ^=1,2,... 

n>0 


are finite. We then expect ps to be positive-definite (at least for some range of 
the parameters) and of trace class. This was proved in [KM] for all the operators 
appearing in Table 2. Moreover, in some cases it is possible to calculate the exact 
integral kernel of the corresponding trace class operator [KM, KMZ, CGuM]. 
An important example in this respect is the family of operators associated to the 
function (63): 

(70) 0,„,„ = e^^ + e^ + e-™"-"^ m,neK>o. 


These operators were called three-term operators in [KM]. In this discussion, m, n 
are positive, real numbers, although in applications to the quantization of mirror 
curves they are often positive integers (like in the quantization of the mirror curves 
to local P(l,n,TO)). Let us first introduce some notation. As in [KM], we will 
denote by <i)b(a;) Faddeev’s quantum dilogarithm [F, FK]. We define as well 


(71) 


g27rarc 

^a,c(a;) = V—7- 

<i)b(a: - i(a -|- c)) 


We then have the following proposition, proved in [KM]. 
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Proposition 3.4. The operator on L^(R) 

Pm,n — Om.n 

is positive-definite and of trace class. Let us define normalized Heisenberg operators 
q, p, satisfying the normalized commutation relation 

(73) [p,q] = (27ri)-^ 


They are related to x, y by the linear canonical transformation, 


(74) 


X = 27rb 


(n -I- l)p -I- nq 
m -I- n -I- 1 ’ 


y = —27rb» 


mp + {m + l)q 


m + n + 1 


so that h is related to b by 
(75) 


h = 


27rb^ 

m + n-\- 1 


In the momentum representation associated to p, the operator pm,n has the integral 
kernel, 


(76) 


Pm.n(p,p') 


^a,c(p) '^a,c{p') 

2b cosh ’ 


In this equation, a, c are given by 


(77) 


mb 

2(m -I- n -I- 1) ’ 


b 

2(to -b n -b 1) 


Once the trace class property has been established for the operators Pm,m it 
can be easily established for operators ps whose inverse Os is a perturbation of a 
three-term operator Om,n by a positive self-adjoint operator [KM]. In addition, it 
can be shown that the operator ps for S' = Fq is also of trace class. This proves the 
trace class property for a large number of operators arising from the quantization 
of mirror curves, including all the operators appearing in Table 2. The trace class 
property of Pm,n and of pfq follows as well from the estimate (68), which was 
established rigorously in these cases in [LSTj. 


3.3. Fredholm determinants from topological strings. In the previous 
sections, we have shown that the operators ps = arising from mirror curves, 
for many toric del Pezzo threefolds, are of trace class. This means in particular 
that their spectral or Fredholm determinant, 

(78) Ss(k,/ i) =det(l-bNps) 

exists and is an entire function of k [Si] . The Fredholm determinant has as a power 
series expansion around k = 0, of the form 

OO 

(79) ^s{^,h) = l+Y,Zs(,N,h)K^, 

where Zs{N^ h) is the fermionic spectral trace^ given by 

(80) Zs{N,h) = Tr{K^{ps)), N = l,2,--- 
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In this expression, the operator A^{ps) is defined by acting on (L^(R)). 
A theorem of Fredholm asserts that, if ps{Pi,Pj) is the kernel of ps, the fermionic 
spectral trace can be computed as a multi-dimensional integral, 

(81) Zs{N,h) = ^ J <iet{ps{pi,Pj)) d^P- 

In turn, the logarithm of the Fredholm determinant can be regarded as a generating 
functional of the spectral traces of ps, since 

(82) Js{K,h) = \og^s{n,h) = - ^ L_^Tr4. 

1=1 

The Fredholm determinant (78) has the infinite product representation [Si] 

OO 

(83) (l-k Ace"'®'*) , 

n—0 

where are the eigenvalues of ps- Therefore, one way of obtaining the spectrum 
of Ps is to look for the zeroes of Sg, which occur at 

(84) K = —e^", n > 0. 

There are very few operators on L^(R) for which the Fredholm determinant 
can be written down explicitly. One family of examples which has been studied in 
some detail are Schrodinger operators with homogeneous potentials V{x) = |a;|®, 
s = 1,2, - • [Vo3, DDT], In this case, the Fredholm determinant is defined as a 
regularized version of 

(85) D { E ) = n (l + f) . 

see for example [Vo2] for a detailed treatment. For these potentials, the function 
D{E) is known to satisfy certain functional equations and it is captured by integral 
equations of the TEA type [DDT]. There is however no closed formula for it. One 
of the surprising (conjectural) results of [GHM] is that the Fredholm determinant of 
the operators ps can be computed exactly and explicitly in terms of the enumerative 
invariants of X encoded in the topological string amplitudes. Moreover, the zeroes 
of the Fredholm determinant are determined by exact quantization conditions. The 
result stated in (31) for the spectrum of /9p2 at S = 27r is a particular example of the 
general conjecture in [GHM]. The resulting quantization conditions are akin to 
those found in conventional Quantum Mechanics (like in, for example, [ZJJ, Vo3]), 
but with an important difference: they involve convergent series, and in particular 
one does not need the apparatus of Borel-Ecalle resummation in order to determine 
the spectrum. 

The conjectural expression of [GHM] for the Fredholm determinant of the 
operator ps requires the two generating functionals of enumerative invariants con¬ 
sidered before, (42) and (52). In order to state the result, we identify the parameter 
K appearing in the Fredholm determinant, with the geometric modulus of X ap¬ 
pearing in (58). We will also write k in terms of the “chemical potential” p, 

( 86 ) 


K = . 
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In addition to the modulus k, we have the “mass parameters” ^j. The flat coordi¬ 
nates for the Kahler moduli space, ti, are related to k and by the mirror map, 
and at leading order, in the large radius limit, we have 

s-l 

(87) - ^aylog^j, f = l, •••,«, 

i=i 

where Ci, aij are constants. As shown in [ACDKV], the WKB approach makes 
it possible to define a quantum mirror map U^h), which in the limit —>■ 0 agrees 
with the conventional mirror map. It is a function of /i, and h and it can be 
computed as an A-period of a quantum corrected version of the differential (45). 
This quantum correction is simply obtained by using the all-orders, perturbative 
WKB method of Dunham. 

We now introduce two different functions of /i. The first one is the WKB grand 
potential, 

F^^(t{h),n) 

h 

+ + A(^, h). 

In this equation, is given by the expression (52). In the second term of 

(88) , the derivative w.r.t. h does not act on the implicit dependence of ti. The 
coefficients bi appearing in the last line are the same ones appearing in (38). The 
function A(^, h) is not known in closed form for arbitrary geometries, although 
detailed conjectures for its form exist in various examples. It is closely related to a 
all-genus resummed version of the constant map contribution Cg appearing in (39). 
The second function is the “worldsheet” grand potential, which can be obtained 
from the generating functional (42), 

(89) b) = F°^ • 

In this formula, B is a constant vector (“B-field”) which depends on the geometry 
under consideration. This vector should satisfy the following requirement: for all 
d, jij and jA such that the refined BPS invariant is non-vanishing, we must 

have 

(90) ^ (_i)B d^ 

For local del Pezzo CY threefolds, the existence of such a vector was established in 
[HMMO]. Note that the effect of this constant vector is to introduce a sign 

(91) 

in the generating functional (42). An important remark is that, in (89), the topolog¬ 
ical string coupling constant gg appearing in (40) and (42) is related to the Planck 
constant appearing in the spectral problem by, 

(92) 9s = ^. 

Therefore, the regime of weak coupling for the topological string coupling constant, 
Ps ^ 1, corresponds to the strong coupling regime of the spectral problem, h ^ 


( 88 ) 


2=1 


jWKB 


27r 


dti 


27r dh 
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1, and conversely, the semiclassical limit of the spectral problem corresponds to 
the strongly coupled topological string. We therefore have a strong-weak coupling 
duality between the spectral problem and the conventional topological string. 

The total grand potential is the sum of these two functions, 

(93) h{pA.K) = + 


and it was first considered in [HMMO]. It has the structure 
(94) 




127rS 


^ ( aijktitjtj- T ^ ( 


i,j,k=l 


f 2TTbi 

[~ir 





where the last term stands for a formal power series in e~*\ whose co¬ 

efficients depend explicitly on h. Note that the trigonometric functions appearing 
in (52) and (42) have double poles when h is a rational multiple of tt. However, 
as shown in [HMMO], the poles cancel in the sum (93). This HMO cancellation 
mechanism was first discovered in [HMO], in a slightly different context, and it 
was first advocated in [KaMa] in the study of quantum curves. 

A natural question is whether the formal power series in (93) converges, at 
least for some values of its arguments. Although we do not have rigorous results on 
this problem, the available evidence suggests that, for real h, J 5 (^,^,h) converges 
in a neighbourhood of the large radius point —)■ oo. However, the series seems to 
be divergent when h is complex. This divergence is inherited from the generating 
functionals (42) and (52). Explicit calculations indicate, for example, that the 
Gopakumar-Vafa generating functional (42) is divergent for complex gs, as first 
noted in [HMMO]. In the one-modulus case, if we write this series as 


(95) F^^{t,gs) = '^ai{gs)e 

e>i 


one finds that, if gs S (D\R, 

(96) log\ai{gs)\ ^ f, £ > 1. 


Therefore, the generating functional (42) is in general ill-defined and, as it stands, 
it can not be used as a non-perturbative definition of the topological string free 
energy^. Finally, note that, after using the explicit expression for the ti in terms of 
pi and the mass parameters, one easily finds that Js(p-, 4, h) is a cubic polynomial 
in /i, plus an infinite series of exponentially small corrections in e”'', At 

large p, the leading behaviour is given by, 

C 

(97) Js(/i,|,h) « —/r>l, 

where C is the constant appearing in (67). 

We can now state the main conjecture of [GHM]. 


Conjecture 3.5. The Fredholm determinant of the operator ps is given by 

(98) 5s(«;, ^> ^) = X! + ^Trin, $.,h)). 


^In some geometries, (42) can be resummed by using instanton calculus, and one obtains a 
convergent function if Qs G C\R [BS]. The resulting function still has a dense set of poles on the 
real line. 
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The sum over n defines a quantum theta function Os (/i, h), 


( 99 ) 




The reason for this name is that, when h = 2?!, the quantum theta function becomes 
a classical, conventional theta function, as we will see in a moment. The vanishing 
locus of the Fredholm determinant, as we have seen in (84), gives the spectrum of 
the trace class operator. Given the form of (99), this is the vanishing locus of the 
quantum theta function. 

It would seem that the above conjecture is very difficult to test, since it gives 
the Fredholm determinant as a formal, inhnite sum. However, it is easy to see 
that the r.h.s. of (98) has a series expansion at large qi in powers of e“^, 

In addition, it leads to an integral representation for the fermionic spectral trace 
which is very useful in practice: if we write Zs{N,h) as a contour integral around 
K = 0, simple manipulations give the expression [HMO, GHM] 


( 100 ) 


Zs{N,h) 



where C is a contour going from e“''^/^oo to e''^/^oo. Note that this is the standard 
contour for the integral representation of the Airy function and it should lead to 
a convergent integral, since Js(/i, S) is given by a cubic polynomial in /i, plus 
exponentially small corrections. Finally, as we will see in a moment, in some cases 
the r.h.s. of (98) can be written in terms of well-defined functions. 

What is the interpretation of the total grand potential that we introduced in 
(93)? The WKB part takes into account the perturbative corrections (in H) to the 
spectral problem defined by the operator ps- In fact, it can be calculated order by 
order in the h expansion, by using standard techniques in Quantum and Statistical 
Mechanics (see for example [MP, H]). The expression (88) provides a resumma¬ 
tion of this expansion at large p. However, the WKB piece is insufficient to solve 
the spectral problem, since as we mentioned above, (88) is divergent for a dense set 
of values of h on the real line. Physically, the additional generating functional (89) 
contains the contribution of complex instantons, which are non-perturbative in h 
[MP, KaMa] and cancel the poles in the all-orders WKB contribution. Surpris¬ 
ingly, the full instanton contribution in this spectral problem is simply encoded in 
the standard topological string partition function. 

The formulae (88), (89), (98) are relatively complicated, in the sense that they 
involve the full generating functionals (52) and (42). There is however an impor¬ 
tant case in which they simplify considerably, namely, when h = 2?!. This was 
called in [GHM] the “maximally supersymmetric case,” since in the closely re¬ 
lated spectral problem of ABJM theory, it occurs when there is enhanced A? = 8 
supersymmetry [GGMj. In this case, as it can be easily seen from the explicit ex¬ 
pressions for the generating functionals, many contributions vanish. For example, 
in the Gopakumar-Vafa generating functional, all terms involving g > 2 are zero. 
A simple calculation shows that the function (93) is given by 
( 101 ) 




1 

871 ^ 


S 




d^Fp 

dtidtj 


47r2 ^ ‘ dU 

i=l 


—Fp{t) + F,{t) + F^^{t), 


where Ko(t)j ^i(t) and Ff^®(t) are the generating functions appearing in (37), (38) 
and (53), with the only difference that one has to include as well the sign (91) in 
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the expansion in e“*'. It can be also seen that, for this value of H, the quantum 
mirror map becomes the classical mirror map, up to a change of sign k —> —k. One 
can now use (101) to compute the quantum theta function appearing in (99). It 
is easy to see that, provided some integrality properties hold for the constant C in 
(67), the quantum theta function becomes a classical theta function. We will now 
illustrate this simplification in the case of local P^. 


Example 3.6. Local has one single Kahler parameter, and no mass param¬ 
eters. The function A(h) appearing in (88) has been conjectured in closed form, 
and it is given by 


( 102 ) 


m 


3Ac{h/Tr) — Ac{3h/Tr) 
4 


where 


(103) Mk) = (l - Q + ^ ^ log(l - 

This function was first introduced in [MP], and determined in integral form in 
[Ha+, HO]. It can be obtained by an appropriate all-genus resummation of the 
constants Cg appearing in (39). In the “maximally supersymmetric case” h = 27r, 
one can write the Fredholm determinant in closed form. The standard topological 
string genus zero free energy is given in (48), (49). The genus one free energies are 
given by [KZ, HKj: 


(104) 


Flit) = I log ^ log (1 + 27z)) , 


Fr(t) = -^log 


1 + 272 


We identify 
(105) 


z = e-^'^ = 


1 


Then, it follows from the conjecture above that 
(^106) 

.:zp 2 (—/v, 27r) = det (1 — Kpp 2 ) 


= exp 1^(271) -h ^ (^Fo(t) - tdtFo(t) + jd^Fo(t)^ + Fi(t) + Ef®(t)| 

where 

(107) ^ ^ (td^Fo(t) - dtFo(t)) , t = ^d^Fo(t), 

and i?2 (z,t) is the Jacobi theta function. The t appearing here is the standard 
modulus of the genus one mirror curve of local P^. In particular, one has that 
Im(r) > 0. It can be easily checked that the function ^ agrees with (29), with the 
identification p. = E. The formulae that we have written down here are slightly 
different from the ones listed earlier in this section, since we are changing the 
sign of K in the Fredholm determinant, but they can be easily derived from them. 

We can now perform a power series expansion around k = 0, by using analytic 
continuation of the topological string free energies to the so-called orbifold pointy 
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which corresponds to the limit z —>■ oo. This is a standard exercise in topological 
string theory (see for example [AGM, ABK]), and one finds [GHM] 

(108) Sp.(«, 2^) = 1 + ^ + 

provided some non-trivial identities are used for theta functions. This predicts 
the values of the very first fermionic spectral traces Zp2(N,2tt), for N = 1,2. 
Interestingly, these traces can be calculated directly in spectral theory, since / 0 p 2 
is the operator pi i defined in (72). The integral kernel of this operator is given in 
(76) (for m = n = 1), and the values appearing in the expansion (108) have been 
verified in this way in [KM, OZ]. 

Finally, it can be easily checked that the function (106) vanishes precisely at 
the values of /r given by the quantization condition (31) (since we are changing the 
sign of K in this formula, the zeros occur at k = e®", n = 0,1, • • •). Note that this 
corresponds to the vanishing locus of the Jacobi theta function appearing in (106). 
In Fig. 4 we show a plot of the Fredholm determinant Sp 2 («;, 27r) as a function of 
K on the real axis. For large, positive values of k, its behavior is dictated by the 
large fi behavior of Jp 2 (/i, 27r), i.e. 

(109) Sp2(K,27r) « exp (log(«;))^| , k > 1. 

On the other hand, for large, negative k, one has the oscillatory behavior, 

(110) Sp2(-K,27r) « expj^ (log(K))^|cos ^TT , k > 1, 

where ^ is again given by (107). □ 



Figure 4. The Fredholm determinant Sp 2 («;, 27r) as a function of 
At, showing the first three zeroes on the negative real axis, corre¬ 
sponding to the first three energy levels —e^", n = 0,1, 2. 

What happens in the general case, when h ^ 27r? One can of course still 
compute the Fredholm determinant from our conjecture (98), but it now contains 
contributions from all the higher genus Gopakumar-Vafa invariants, and from the 
all-orders NS free energy. In order to study its properties, one possibility is to 
expand it in power series around the “semiclassical” case h = 27r. This was done in 
[Gr], and the resulting expansion puts on a rigorous footing the 1/A expansion of 
the “non-perturbative partition function” studied in [EM]. One can also use (98) 
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to determine the spectrum by looking at the vanishing locus of the Fredholm deter¬ 
minant for general h. By using the large radius expansion of the grand potential, 
one can write an exact quantization condition defined by an explicit power series. 
As shown in [GHM] in many examples, the predictions for the spectrum obtained 
from this exact quantization condition are in perfect agreement with direct numer¬ 
ical calculations (as done for example in [HW]). It was found in [WZH] that, 
when the mirror curve has genus one, the quantization condition of [GHM] can 
be written down in closed form in terms solely of the NS free energy. The results 
of [WZH, GHM], when put together, imply that the generating functional (52), 
capturing the NS limit of the topological string, is closely related to the generat¬ 
ing functional (42), capturing the standard topological string. This connection has 
been further explored in [SWH] and finally established in [GG] in many cases. 

For a generic value of h, it is difficult to extract analytic results for the spectral 
traces from (4), since the analytic continuation of the generating functionals (42), 
(52) to the orbifold point for arbitrary h is not known. However, one can still obtain 
very precise numerical results. The basis for this is the formula (100). By expanding 
8.t large /i, the fermionic spectral trace Zs{N,h) can be evaluated as an 
infinite sum of Airy functions and their derivatives. Let us spell this out in some 
detail for local P^. One can write. 


( 111 ) 


eJp2(MA) = ^ai,„/r”e-''^, 

Z,n 


where 

(112) = + + 


The values of C{ti), B{h) can be read off from the general expression (93) and are 
given by [GHM] 


(113) 



B{h) = 


TT 

2h 


n 

IOtt 


Then, assuming good convergence properties for the expansion of Jp 2 (/x, h), we can 
exchange it with integration in (100). We obtain in this way. 


(114) 


Z-p2 (A^, A) 


e^(^) 

(C(^ 


l,n 



f N + l-B{h) \ 

V J 


where Ai(z) is the Airy function. Note that n takes non-negative integer values, 
but I is of the form 3p -I- Girq/H, where p, q are non-negative integers. The leading 
behavior of (114) in the limit N ^ oo, h fixed, is given by the first term in the 
r.h.s., which is an Airy function 

(115) JV»1. 

The additional terms in the r.h.s. of (114) give an infinite series of exponentially 
suppressed corrections to (115). Remarkably, this series seems to be convergent, and 
it should be regarded as the analogue in this theory of the Rademacher expansion in 
number theory. It also produces highly accurate numerical answers for the spectral 
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traces Z-p 2 (N,h). Using this procedure, one obtains for example, 
(116) 


Zp2 


0.4604521481728325977904889856168747087632124207... 


This can be compared to the analytic result obtained by integrating (76), 

(117) Trpi,i = I) = 1 exp , V = 2Im [u, (e^^/^)) . 


The agreement is quite remarkable, and the precision increases with the number 
of terms retained in the expansion of the grand potential. Many tests of the main 
conjecture (98), for various local del Pezzo CYs, can be found in [GKMR]. In all 
cases, one finds that (98) captures in exquisite detail the spectrum and the spectral 
traces of the corresponding operator. 

It follows from (115) that, at large N and h fixed, the spectral traces display 
the scaling 

(118) \ogZp2{N,H)^-^N^/^H^/^, iV>l. 

9 

This can be also derived as a consequence of the leading WKB scaling (24), see 
for example [MP] for a detailed explanation of this type of relation. The above 
behaviors, (115) and (118), are not restricted to local P^, but hold for all local del 
Pezzo CYs: the leading term in the large N expansion of the fermionic spectral 
trace is always an Airy function of the form (115), which leads to the scaling. 
These behaviors are typical of a theory of N M2 branes [KT], and they are very 
similar to what happens in ABJM theory [DMP, FHM, MP]. A universal Airy 
behavior for topological strings has been recently discovered in [AYZ, CS], and 
it would be interesting to understand its precise relation to what is found in this 
formalism. 


4. Prom spectral theory to topological strings 


4.1. Random matrices from operators and topological strings. An im¬ 
portant consequence of the conjecture of [GHM] is the following. As we have seen, 
the Fredholm determinant involves both the WKB grand potential, which is cap¬ 
tured by the NS free energy, and the worldsheet instanton grand potential, which 
is essentially given by the topological string free energy. It turns out that, due to 
their functional form, there are appropriate scaling limits in which only one of these 
ingredients contributes. We will focus on the’t Hooft-like limit 


(119) 


h —>■ 00 , H —)■ oo. 


^ ( fixed. 

n 


For simplicity, we will also assume that the mass parameters scale in such a way 
that 


( 120 ) 


nij = i 


2-t: jh 


j = 1,... ,S- 1, 


are fixed (although other limits are possible, see [MZ] for a detailed discussion). 
It is easy to see that, in this limit, the quantum mirror map becomes trivial, and 
the approximation (87) is exact. The grand potential has then the asymptotic 
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expansion, 


(121) 


where 

Jo (C,»T^) 

(122) 

Jf (C,m) 

Jg (C,"i) 

In these equations. 

(123) 



3=0 


Fq (t) + 471^ ^ 6f + 147r"‘ylo (m) 


IGtt^ 


(C, m) = ylg (m) + (47T^ f^ 2 - Cg) , g > 2. 


ti = 27rciC - X! 

i=i 


and Fg (t) are the standard topological string free energies as a function of the 
standard Kahler parameters t, after turning on the B-field as in (91). We have also 
assumed that the function A {^,h) has the expansion 

OO 

(124) A{tn) = J2Mm)n^~"^, 

3=0 

and this assumption can be tested in examples. 

Therefore, in the ’t Hooft limit (119), the grand potential has an asymptotic 
expansion which is essentially equivalent to the genus expansion of the standard 
topological string. According to our conjecture, this expansion gives a concrete 
prediction for the ’t Hooft expansion of the fermionic spectral traces, thanks to 
the equality (100). Indeed, let us consider the corresponding ’t Hooft limit of the 
fermionic spectral traces, 

N 

(125) h —>■ OO, N —7> OO, — = A fixed. 

h 

Then, it follows from (100) that one has the following asymptotic expansion, 


(126) log Zs{N,fi) = 

3>0 

where -FgiX) can be obtained by evaluating the integral in the right hand side of 
(100) in the saddle-point approximation and using just the ’t Hooft expansion of 
Js(Mj^)^)- At leading order in h, this is just a Legendre transform, and one finds 

(127) J-o(A)=J^(C,m)-AC, 


evaluated at the saddle-point given by 


(128) 


dC ■ 


The higher order corrections can be computed systematically. In fact, the formalism 
of [ABK] gives a nice geometric description of the integral transform (100) in the 
saddle-point approximation: the functions ^g{\) are simply the free energies of the 
topological string on X, but on a different frame. This frame corresponds to the 
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so-called conifold point of the geometry. In particular, A turns out to be a vanishing 
period at the conifold point. 

Therefore, the conjecture (98), in its form (100), provides a precise prediction 
for the’t Hooft limit of the fermionic spectral traces: they are encoded in the stan¬ 
dard topological string free energy, but evaluated at the conifold frame. It turns 
out that this prediction can be tested in detail. The reason is that, at least in some 
cases, the fermionic spectral traces can be expressed as random matrix integrals, 
and their ’t Hooft limit can be studied by using various techniques developed for 
matrix models at large N. In the case of the three-term operators (70), the ex¬ 
plicit expression for the kernel appearing in (76) can be combined with Fredholm’s 
theorem (81) to obtain the following expression for their fermionic spectral traces 
[MZ], 

(129) 


Zm,u{N, h) 


1 

m 




^ n.<,4sinh(^)' 

Hij 2 cosh + iTrC'™,^^ 


This can be regarded as an 0(2) matrix model [K, EK], with a potential depending 
on h and given by 


(130) 


V (u, h) 


-h\og 




In [MZ] it is shown in detail that this potential admits an asymptotic expansion at 
large S, which can be used to obtain the’t Hooft expansion of Zjn,n(N, h). In par¬ 
ticular, when m = n = 1, one finds a description of the all-genus topological string 
free energy of local as the asymptotic ’t Hooft expansion of this matrix integral. 
Detailed calculations show that the expansion obtained in this way agrees exactly 
with the predictions of (100). In fact, the conjecture “explains” many aspects of 
topological string theory at the conifold point, for many toric CY geometries. For 
example, it has been known for some time that the leading terms in the expansion 
of the conifold free energies are universal [GhV], and that they agree with the 
all-genus free energy of the Gaussian matrix model. Since matrix models like (129) 
are deformations of the Gaussian one, this observation is a consequence of our con¬ 
jecture (or, equivalently, this observation must hold in order for the conjecture to 
be true). 

Note that, as shown in Fig. 5, different regions of the moduli space of a toric 
CY are associated to different expansions of the spectral quantities: the large radius 
region defines the grand potential at large The orbifold point is the appropriate 
one to calculate the Taylor expansion of the Fredholm determinant at the origin, 
as in (108), while the conifold point encodes the ’t Hooft limit of the fermionic 
spectral traces. 


4.2. Non-perturbative topological strings. As we mentioned above, the 
total free energy of the topological string (40) is hopelessly divergent. This is also 
the case for the total free energy at the conifold frame. Let us consider the formal 
power series in H appearing in the r.h.s. of (126). As it happens at large radius, 
the functions J^g(X) have a common, finite radius of convergence. For each A inside 
this radius of convergence, the resulting series in h diverges factorially, like 

(131) J-g(A) ~ (25)!(A(A))-2^ 5»1- 
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large radius 


K) 



Figure 5. Different points in the moduli space lead to different 
expansions for the spectral quantities. 

Therefore, the expansion (126) is asymptotic and does not define a function of A, 
h. However, we have now found a quantity, the fermionic spectral trace, with the 
following properties: 

(1) It is rigorously well-defined. This is a consequence of the trace class prop¬ 
erty of PS- 

(2) Its asymptotic expansion, in the ’t Hooft limit, reproduces exactly the 
genus expansion of the topological string in the conifold frame. 


mirror curve 


trace class operator 

W{x,y) 

quantization 

P 


fermionic 

trace 


topological string 


matrix model 

log Z{N, h)^Yu 

9>0 

’t Hooft limit 

Z{N, h) 


Figure 6. Given a toric CY threefold, the quantization of its mir¬ 
ror curve leads to a trace class operator p. The standard topological 
string free energy is obtained in the’t Hooft limit of its fermionic 
traces Z{N,h). 

It follows that the fermionic spectral traces provide a non-perturbative comple¬ 
tion of the topological string, in the conifold frame. The results of [GHM, MZ, 
KMZ] lead then to an elegant approach to non-perturbative topological strings: 
given a toric CY manifold X, one considers the quantization of its mirror curve, 
which leads to a trace class operator. The standard, all-genus topological string 
free energy emerges as an asymptotic expansion of the fermionic spectral traces of 
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the operator, in a ’t Hooft-like limit. This procedure is sketched in Fig. 6. Note 
that, from the point of view of the spectral problem, the conventional topological 
string emerges at strong coupling 3> 1. 

Non-perturbative completions are not unique. One way to obtain well-defined 
answers from asymptotic series is to use Borel resummation. In favorable cases, 
when the series is Borel summable, this procedure is unambiguous and leads to 
a function whose asymptotic expansion reproduces the original asymptotic series 
(see for example [M2] for a review). In the case of topological string theory in the 
conifold frame, it turns out that, in many examples (like local P^), the genus expan¬ 
sion is Borel summable for real h and real A. The result of the Borel resummation 
however turns out to be different from the non-perturbative completion defined by 
the fermionic spectral traces. For example, in the case of local P^, one finds that 
[CMS] 

(132) log 27r) = -2.197217... 

to be compared with the value obtained from spectral theory (see (108)), 

(133) log Zp 2 (1, 27r) = - log(9) = -2.197224... 

The difference between the Borel resummation and the answer from spectral the¬ 
ory is due to exponentially small, non-perturbative effects associated to complex 
instantons. In the context of the theory of resurgence, this indicates that the per¬ 
turbative genus expansion has to be generalized to a trans-series including the 
non-perturbative effects in the form of additional asymptotic series. Such trans¬ 
series can be produced by using the powerful techniques of [CESV], which are 
based on a generalization of the holomorphic anomaly equation of [BCOV2] to 
the non-perturbative realm. One can now wonder whether the non-perturbative 
information contained in the fermionic spectral traces can be recovered from the 
Borel-Ecalle resummation of the full trans-series. In [CMS], strong numerical 
evidence has been given that this is the case relating in this way the “exact” non- 
perturbative approach of [GHM] and the resurgent analysis of non-perturbative 
effects in [CESV]. 

The fermionic spectral trace Z{N,h) provides a non-perturbative completion 
of the topological string in the conifold frame, and it actually promotes it to an 
entire function of N, as first noted in [CGM]. In the large radius frame, we can 
also promote the topological string partition function to a entire function on moduli 
space, which is nothing but the spectral determinant itself. More precisely, in the 
’t Hooft regime (119), we have the asymptotic behavior 

(134) logSs(/c,|,S)^J^Hooft(t)^ 

In this equation, the ti are related to m as spelled out in (123). The r.h.s. is 
given in (121) and (122), and it involves the genus expansion of the topological 
string in the large radius regime. Therefore, we can promote the total topological 
string partition function in the large radius frame, to an entire function of moduli 
space. It was pointed out by Witten in [W] that this partition function should be 
thought as a wavefunction on the moduli space of the CY. In this view, the genus 
expansion should be regarded as a sort of WKB approximation to this wavefunction. 

^Non-perturbative effects are explicitly included in the NS contribution to the grand potential, 
but this contribution does not have the standard form of a trans-series in the string coupling 
constant. 
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However, the WKB solution has branch cuts and singularities which are an artifact 
of the approximation, and are not present in the exact answer, which is sometimes 
an entire function. In [MMSS], such a phenomenon has been argued to occur as 
well in string theory, where perturbative wavefunctions should be “smoothed out” 
by non-perturbative effects. It is tempting to interpret (134) in this light, and to 
regard the spectral determinant as the exact, entire wavefunction on moduli space 
postulated in [W]. As in the WKB approximation, the branch cuts and conifold 
singularities appearing in the perturbative free energies of the topological strings 
should be regarded as artifacts of the genus expansion, and they are smoothed out 
in the final answer. 


4.3. Number-theoretic aspects. It has been observed in [RV, DK] (see 
also [MOY]) that the values of the Kahler parameters of many toric CYs at the 
conifold point are given by special values of the dilogarithm function. It turns out 
that these number-theoretic properties of the CY periods at the conifold point are 
consequences of the main conjecture (98). According to this conjecture (98), the 
topological string amplitudes are encoded in the spectral properties of the operator 
obtained by quantizing the curve. On the other hand, as shown in [KM], this 
operator is closely related to the quantum dilogarithm. When these two statements 
are put together, the results of [RV, DK] follow. 

To see this in more detail, note that (127) implies that 

(135) ^ ^ -C, 

It follows from (126) that the l.h.s. of this equation can be calculated by studying 
the large N limit of the fermionic spectral traces. For the operators Om,n, it can 
be computed from the matrix integral (129), and the result is [MZ] 


(136) 

In this equation. 


dX 


A=0 


m -\- n + 1 
271 ^ 


D{-q 


m+1 


Xn 


(137) 


q = exp 


m n 1 


Xk 



and D{z) is the Bloch-Wigner function defined by, 


(138) D{z) = ImLi 2 ( 2 :) + arg(l - z) log \z\, 


where arg denotes the branch of the argument between — tt and tt. Now, the point 
A = 0 is the conifold point, and since C is related to the Kahler parameters ti by 
(123), one can evaluate these parameters at the conifold point in terms of (136). 

The simplest example of this situation is again local P^. In this case there is a 
single Kahler parameter t, related to C by 

(139) t = 6 ttC. 

The conifold point occurs at z = —XjTl, which is at the boundary of convergence 
of the series in (28). The value of t at this point, which we will denote by t^, can 
be in principle obtained by evaluating the mirror map (47), 

(140) t, = -log(l/27)-Wi(-l/27). 
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Here, we are dropping an imaginary piece ±7ri coming from the log. According to 
(136), (140) can be evaluated in terms of the Bloch-Wigner function as 



(141) 


This is clearly a non-trivial number-theoretic property of the period. In the frame¬ 
work of [GHM, MZ], it follows ultimately from the fact that the operator / 0 p 2 
involves the quantum dilogarithm [KM]. The result (141) turns out to be true, 
and it has been proved in [RV, DKj. Similar identities can be obtained for other 
CYs, and some of them seem to be new (see [CGM2, CGuM] for examples in¬ 
volving mirror curves of genus two.) 

The fact that the spectral theory realization of the topological string “explains” 
some number-theoretic properties of the periods seems to indicate that this non- 
perturbative completion is capturing essential aspects of the theory. 


5. Outlook 


In this paper we have reviewed the correspondence between spectral theory 
and topological strings presented in [GHM]. This correspondence leads to a new 
family of exactly solvable problems in operator theory, providing explicit expressions 
for the corresponding Fredholm determinants. At the same time, it leads to new 
insights on topological strings on toric CY threefolds, and to a non-perturbative 
definition of these string theories in the spirit of large N string/gauge dualities. 

In this review, we have focused on mirror curves of genus one. The case of genus 
zero is certainly special, since the corresponding operators do not have discrete 
spectra. However, the ideas of [HMMO, GHM] can be extended to this case in a 
natural way, as it has been noted recently in [H, H02, KrMk]. More interesting, 
and more difficult, is the extension to mirror curves of higher genus. This has been 
done in [CGM2] and further explored in [CGuM]. When the mirror curve has 
genus gs, its quantization leads to gs different trace class operators. It is possible to 
define a generalized Fredholm determinant depending on g-^ moduli. There are also 
natural generalizations of the fermionic spectral traces, which now depend on g-^ 
different non-negative integers Ni, i = I,-- - ,g-s, and provide a non-perturbative 
definition of the topological string in a ’t Hooft-like limit. However, in contrast 
to the situation for quantum integrable systems, there is a single quantization 
condition, which determines a codimension one locus in moduli space through the 
vanishing of the Fredholm determinant 

The results presented here have important consequences for the study of a class 
of integrable systems, called cluster integrable systems, which can be associated to 
local CY manifolds. These systems were constructed by Goncharov and Kenyon in 
[GoKe] and include in particular the relativistic Toda lattice [EFS, Mar, FM]. 
The number of Hamiltonians in these systems is equal to the genus of the mirror 
curve. In the case of curves of genus one, the operator Os is nothing but the 
(single) quantum Hamiltonian of the Goncharov-Kenyon system. Therefore, the 
results of [GHM] provide a conjectural solution to the quantum integrable system 
in the genus one case. By using the formulation of the quantization condition of 
[GHM] presented in [WZH], an exact quantization condition for the spectrum of 
the general Goncharov-Kenyon system has been conjectured in [HM, FrHaMa]. 
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There are clearly many avenues to explore in the future. Although we have 
presented some detailed results for the spectrum of the operators, one should ad¬ 
dress the construction of the eigenfunctions, which is closely related to the open 
string sector. A first step in this direction has been made in [MZ2], where a gen¬ 
eral proposal has been made in the case of genus one curves, and an exact solution 
has been written down explicitly in the maximally supersymmetric case of local 

X P^. However, much more work is needed in order to have a general solu¬ 
tion for the eigenfunctions. Eventually, one should provide a rigorous proof (or at 
least a derivation at the physics level of rigor) of the quantization condition. This 
looks difficult, since, for this type of operators, we do not even have a semi-rigorous 
method to calculate instanton corrections. Developing these tools is an important 
task for the future. 

On the topological string side, one limitation of the framework we have pre¬ 
sented is that it does not accommodate the general refined theory, but only the two 
one-parameter specializations corresponding to the standard and the NS limits. 
However, one lesson we have learned from the developments reviewed in this paper 
is that these two one-parameter cases are not independent: in order to solve the 
spectral problem, we need both generating functions, (42) and (52), and each one 
of them can be regarded as a non-perturbative correction to the other. Maybe we 
should think about the refined topological string with general (£ 1 , 62 ) parameters 
as a non-perturbative object, similar to the {p,q) string, instead of looking for a 
perturbative worldsheet description. 
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